Fee breathing instability in BaBiC>3 from first principles 
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We present first-principles density-functional calculations using the local density approximation to 
investigate the structural instability of cubic perovskite BaBi03 . This material might exhibit charge 
disproportionation and some evidence thereof has been linked to the appearance of an additional, 
fourth peak in the experimental IR spectrum. However, our results suggest that the origin of this 
additional peak can be understood within the picture of a simple structural instability. While the 
true instability consists of an oxygen-octahedra breathing distortion and a small octahedra rotation, 
we find that the breathing alone in a fec-type cell doubling is sufficient to explain the fourth peak 
in the IR spectrum. Our results show that the oscillator strength of this particular mode is of the 
same order of magnitude as the other three modes, in agreement with experiment. 

PACS numbers: 77.80.-e, 71.20.-b, 63.20.-e, 78.30.-j 



Over the last few decades, BaBi03 has been a subject 
of continuing interest for several reasons. First, this ma- 
terial is an end member of the series (Ba,K)(Pb,Bi)C>3, 
which includes one of the few high-temperature super- 
conductors that do not contain copper. Second, the for- 
mal valence of Bi for the cubic perovskite BaBiC>3 is 4+, 
suggesting that the Bi should exhibit charge dispropor- 
tionation 2Bi 4+ — > Bi 3+ + Bi 5+ leading to a charge- 
ordered state. A better understanding of such ordered 
states could lead to the design of new charge-ordered fer- 
roelectrics, a class of materials that has recently attracted 
considerable attention. 1 Experimental evidence for static 
charge disproportionation in BaBiC>3 has been linked to 
the interpretation of an additional, fourth strong peak in 
its IR spectrum^ai^ while group theory predicts only 
three peaks for the cubic perovskite structure. 

The goal of this paper is to show that not only 
charge disproportionation, but also a much simpler 
mechanism — i.e. a structural breathing instability — can 
account for all four peaks in the experimental IR spec- 
trum of BaBi03. To this end, we first perform density- 
functional theory (DFT) calculations to investigate the 
structural instability of cubic perovskite BaBiC>3. There- 
after, we calculate the IR active phonons and the Born 
effective charges, which in turn allows us to determine 
the IR oscillator strength of all four modes mentioned 
above. We would also like to point out that for our cal- 
culations we assumed a simple fee cell-doubling breathing 
distortion, which is sufficient to warrant non-zero oscil- 
lator strengths for all four modes. The full structural 
instability, however, consists not only of a breathing of 
the oxygen octahedra, but also of a small nearly-rigid 
rotation thereof AL£ 

Our first-principles density-functional-theory calcula- 
tions were performed using the VASP package 9 and PAW 
potentials^ with a 520 eV kinetic-energy cutoff. Note 
that we included the Ba 5s 2 and Bi 5c£ 10 states as valence 
states. The exchange-correlation energy functional was 
evaluated within the local density approximation (LDA) 
as parameterized by Perdew and Zunger^l Brillouin-zone 



FIG. 1: Distortion of the oxygen octahedra. Small circles de- 
note oxygen atoms and big circles denote bismuth atoms, re- 
spectively. In three dimensions, this breathing distortion pro- 
duces a checkerboard pattern of adjacent smaller and larger 
octahedra, doubling the primitive simple cubic unit cell to fee. 



integrals were approximated with the Monkhorst-Pack 
scheme on a 14x14x14 mesh of special k-points. 12 Spin- 
orbit effects were not taken into account, since previous 
calculations showed that their effect on the band struc- 
ture around the Fermi energy is small^ 

The high-temperature phase of BaBi03 above 750 K 
has the simple cubic perovskite structure with five atoms 
per unit cell£ The simple cubic lattice constant of 
4.328 A for our calculations was derived from reported 
experimental values. In this structure, BaBiOs has an 
odd number of electrons per unit cell and therefore is 
metallic. However, lower-temperature experiments and 
calculations show that the unit cell for the low temper- 
ature phase is monoclinicj§iL& and that the material is 
a semiconductor with an estimated gap of 0.2 eVii^ The 
monoclinic phase results from a distortion of the cubic 
cell consisting mainly of a breathing of the oxygen octa- 
hedra and a small nearly-rigid rotation thereof. As men- 
tioned above, for our investigation we limit ourselves to 
the pure breathing of the octahedra, which is enough to 
explain the appearance of all four peaks in the IR spec- 
trum. In the breathing distortion, the oxygen atoms can 
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move along the Bi-0 bond to produce two inequivalent 
Bi sites with different Bi-0 bond lengths, without break- 
ing any point symmetries. This results in a doubling of 
the unit cell to ten atoms, following the three dimensional 
checkerboard pattern sketched in Fig. ^ The resulting 
fee lattice has a conventional lattice constant of 8.655 A. 
This distortion has previously been proposed by Uchida 
and coworkers. 2 

Figure [2] depicts the resulting double-well potential 
when the oxygen atoms are moved along the Bi-0 bonds 
in the vicinity of the high symmetry oxygen position of 
1/4, marked as A. The different curves correspond to dif- 
ferent dense k-meshes. Our results show that a 8 x 8 x 8 
mesh, i.e. 60 k-points in the irreducible part of the Bril- 
louin zone (IBZ), which is usually sufficient for calcula- 
tions on similar perovskites, does not provide converged 
results. A very dense mesh of 14 x 14 x 14 k-points (i.e. 
280 k-points in IBZ) is necessary to describe this sen- 
sitive breathing distortion correctly. This is partly due 
to the fact that — as we will see below — the system is 
close to the metal/semiconductor transition. The double 
well has minima of depth 6.7 meV per 10-atom unit-cell 
at a displacement of ±0.00507 (B), corresponding to the 
symmetry-equivalent choices of breathing a given octahe- 
dron either in or out. The shape of our double- well po- 
tential agrees well with LMTO calculations for breathing 
where a rotation of 13° has been frozen in>i£ Other calcu- 
lations find for the double-well without rotation a depth 
of 20 meV at a larger displacement of 0.07 A (LMTO)£ 
and a depth of 50 meV at 0.06 A (LAPW) Lichtenstein 
and coworkers, including combined tilting and breathing, 
find a minimum of 14 meV at 0.055 A with 8.5° rotation^ 
Note that the experimental values at 150 K are a rotation 
of 11.2° and a displacement of 0.085 A. 

Henceforth, we will abbreviate in this paper our 
calculated optimal oxygen displacement as ±0.00507« 
±0.005, keeping in mind that all calculations were actu- 
ally done at the exact minimum position. At the min- 
ima, the oxygen atoms are displaced 0.04A. As a result, 
the simple cubic Bi-0 bond length of 2.164 A splits into 
2.120 A and 2.208 A for Bi(l)-0 and Bi(2)-0. We also 
studied the shape of the double-well potential as a func- 
tion of pressure. It turns out that pressure decreases the 
depth of the double well and moves its minimum closer 
to the high-symmetry position, thereby stabilizing the 
structure. 

In order to investigate the effect of the breathing upon 
the gap, we calculated the band structure and DOS 
for different oxygen displacements as depicted in Fig. 
It can be seen that at the minima of the double-well 
(B) a small band overlap along W-L is present, rather 
than a gap. This finding is in agreement with previous 
calculations 13 and the discrepancy is consistent with the 
well-known tendency of LDA to underestimate the gap. 
For the larger displacement according to point C, the 
gap is clearly open. It can also be seen that the degener- 
ate band along W-L in the simple cubic structure splits 
as the oxygens move. The splitting is almost constant 
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FIG. 2: K-mesh dependence of the total energy of fee BaBiC>3 
as a function of oxygen displacement from the high symmetry 
position of 1/4 (point A) given in relative Cartesian coordi- 
nates. Results are plotted for meshes from 8x8x8 (i.e. 60 
k-points in the IBZ) up to 16 x 16 x 16 (i.e. 408 k-points in the 
IBZ). Convergence is reached with a 14 x 14 x 14 mesh, which 
produces results indistinguishable from any denser mesh. The 
minima of the double- well potential are at ±0.005 (B), result- 
ing in octahedra of two different sizes. 



along the band and of the order of 1 eV. Calculations 
by Mattheiss and Hamann, including the breathing and 
tilting of the octahedra, show that this splitting is pri- 
marily due to the breathing^ and is therefore captured 
correctly within our approach. In addition, the bands at 
K repel each other. Both effects together provide a mech- 
anism to open the gap in agreement with the observed 
semiconducting character of BaBiOa. 14 

We next used a frozen-phonon approach to calculate 
the zone-center phonon frequencies and eigenvectors for 
the minimum-energy fee structure with oxygen displace- 
ments of ±0.005, corresponding to point B in Fig. [3 The 
optical modes in this case are 



r ^ = A lg + E g + T 2u + 2T 2g + AT lu + T lg 



(1) 



four of which are IR active (T\ u ). The symmetry adapted 
calculation gives for these four modes the frequencies 
Wl = 444 (438), lu 2 = 198 (231), lu 3 = 157 (136), and 
L04 = 32 (98) cm -1 , where the experimental values in 
parenthesis are taken from Ref. |j. Considering that 
we neglected the rotational distortion, the agreement of 
the high-frequency modes is reasonable, while the low- 
frequency mode at 31 cm -1 differs considerably from the 
experimental value. 

Next, we calculated the Born effective charge tensors 
%k up usm g a Berry-phase schemeiiii^ This calculation 
could not be carried out directly for the minimum-energy 
fee structure, as it is metallic due to a small band over- 
lap, as described above. We computed the Born effec- 
tive charges for the semiconducting structure obtained 
by displacing the oxygens by an additional 0.03A to the 
experimental displacement of 0.08A, since at this point 
the gap has just opened. This procedure is justified by 
the fact that the phonon corresponding to the oxygen 
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— 0.000(A) 

— 0.005(B) 

— 0.010(C) 




When the unit cell is doubled as described above, the 
zone-boundary i?-point phonons arc folded back to the 
zone center. This leads to a mixing of one of the folded- 
back phonons (with alternating Bi displacements relative 
to their surrounding oxygen octahedra) and the IR active 
zone-center phonons. In the five-atom cell this phonon 
has no oscillator strength since the contributions of neigh- 
boring Bi cancel. In the doubled cell, the octahedra of 
different sizes create different environments for the bis- 
muth atoms sitting in their centers and the contributions 
no longer cancel. As a result, the oscillator strength of 
this particular mode becomes non-zero, which explains 
the additional, fourth peak in the experimental IR spec- 
trum. Note that materials such as BaTi03 or SrTiO.4, 
that have a single valence for the Ti ion (4+), only show 
three strong peaks in the IR spectrum — as expected from 
group theory. 

As the last step, we calculate the oscillator strengths of 
the IR active modes. The oscillator strength Si of mode 
i can be calculated as 



FIG. 3: Band structure and DOS of BaBiOs for the oxy- 
gen displacements A, B, and C of Fig. [5] Calculation A 
corresponds to the simple cubic perovskite cell and reflects 
the metallic character expected. The oxygen displacement of 
0.005 (B) splits the degenerate band along W-L and decreases 
the DOS around the Fermi energy. The results are plotted 
with respect to the Fermi energy and the DOS is given in 
states/eV. 



TABLE I: Born effective charges Z% a a for atom k in units of 
|e| for fee BaBiOs. All Z| a g are diagonal. Bi(l) and Bi(2) 
are located in the center of a larger and smaller octahedra, 
respectively. For the O-atoms we refer to ZTi and Z*[ , denoting 
the components parallel and perpendicular to the Bi-O bond. 



Atom 


"^fc.lli -^fc,22i ^k,33 


Atom 






Bi(l) 


+4.78 


0(1-6) 


-4.55 


-1.85 


Bi(2) 


+6.22 








Ba(l), Ba(2) 


+2.75 









displacement is non-polar and the Born effective charges 
are expected to be smooth near the metal/semiconductor 
transition. Indeed, we investigated the dependence of 
the effective charges upon the oxygen displacement by 
performing calculations for a larger displacement corre- 
sponding to 0.010 (point C in Fig. El- Our results show 
that the effective charges for the displacement 0.009 and 
0.010 differ only by approximately 1%. All Z% a/3 are 
found to be diagonal and the results are collected in Ta- 
ble Q] The values for the effective charges are reminiscent 
of similar perovskites, such as BaTi03ii^ The splitting 
of Z* for Bi(l) and Bi(2) is quite striking, however, note 
that this "Born effective charge disproportionation" is 
not necessarily the result of static charge transfer. 

For the simple cubic perovskite structure group the- 
ory predicts three 3-fold degenerate IR active phonons. 2 



a 1 p,k 



7* 



6a (M) 



-1 2 



(2) 



where £(k,i) is the displacement of atom k according 
to the eigenvector of the ith phonon mode, ink is the 
mass of atom k, and Z£ a o are the Born effective charges. 
Our results for the oscillator strengths are S\ = 2.097, 
S 2 = 0.457, S3 = 0.661, and S 4 = 0.617. The three 
modes Si (i — 1,2,4) have been associated with the sim- 
ple cubic perovskite structure and can be identified in 
the experimental spectrum£2ii£ The oscillator strength 
of the fourth mode S3 — 0.661 is of the same order of 
magnitude as the ones associated with the simple cu- 
bic perovskite structure. Indeed, as a result, a fourth 
peak is clearly visible in experiments. The fact that our 
calculations predict an oscillator strength for this mode 
comparable to the other modes goes beyond pure group 
theory, which cannot make any predictions about the os- 
cillator strength. As our main result, we conclude that 
DFT calculations for a pure fec-type breathing can ex- 
plain the existence of all four peaks in the IR spectrum 
of BaBi0 3 . 

Further analysis of 53 reveals that its value is al- 
most independent of the Born-effective-charge splitting 
between the two inequivalent Bi atoms. If we remove 
the splitting between — 4.78 and = 6.22 by 

hand and replace these values with the average of 5.5, 
the oscillator strength changes only by a few percent. It 
follows that the oscillator strength is mainly determined 
by the eigenvector of the corresponding phonon mode. 
Furthermore, we find that about half of the oscillator 
strength is a result of the asymmetry in displacement of 
the Bi(l) and Bi(2) atoms. The rest is mostly due to 
the displacement of the four oxygen atoms in the plane 
perpendicular to the Bi movement. 

To conclude, we have performed first-principles DFT 
calculations to investigate the electronic structure and 
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structural instability of cubic perovskite BaBiC>3. We 
investigate the effect of a fcc-like unit-cell doubling pro- 
duced by a breathing distortion of the oxygen octahedra. 
The oxygen atoms are thereby displaced 0.04A along the 
Bi-0 bond, which generates smaller and larger octahe- 
dra throughout the crystal. The distortion is very small, 
and yet we find that properties such as the band gap or 
the Born effective charges are very sensitive to it. This 
cell-doubling distortion turns-on the oscillator strength 
of an additional fourth mode, observed as an additional 
strong peak in experiment. In particular, our calculations 
show that the oscillator strength of this certain mode is 



of the same order of magnitude as the strengths of the 
remaining three modes. In future work, we would like to 
improve our results by including not only the breathing 
distortion, but also the rotation of the octahedra. Fur- 
thermore, although BaBiC>3 is not considered a strongly 
correlated system, it might be interesting to investigate 
the effect of LDA+U in conjunction with the structural 
instability. 
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